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Анал1з методу усунення «витоку» спектра 
при длагностиц! систем турбонаддуву судових дизелв 


Дан анализ численного метода устранения эффекта «утечки» спектра, полученного в результате 
дискретного преобразования Фурье (ДПФ). Практическое применение метода заключается в том, что он 
позволяет оценивать основные параметры измеряемого периодического сигнала (амплитуды И частоты) по 
его спектру. Анализируемый в статье метод является альтернативой применения оконных функций, 
подавляющих, но не исключающих эффект «утечки» спектральной мощности. В отличие от применения 
ОКОННЫХ функций, представленный в статье метод позволяет практически полностью восстановить 
истинные значения амплитуды и частоты сигнала по данным дискретного спектра. В статье анализируется 
работа предлагаемого численного метода устранения эффекта «утечки» спектральной мощности при 
различной степени зашумленности исходного периодического сигнала. В качестве примера взят 
анализатор спектра виброакустического сигнала турбины наддува судового дизеля, построенный на базе 
контроллера АТхтега256АЗ. 

Ключевые слова: спектральный анализ, эффект «утечки» ДПФ, судовой дизель, 
турбонаддув, диагностика. 
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Наведений аналз чисельного методу усунення ефекту «витоку» спектра, отриманого в результати 
дискретного перетворення Фур’е (ДПФ). Практичне застосування методу полягае в тому, що з його 
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допомогою стае можливою оцнка основних параметрв вимрюваного перюдичного сигналу (ампттуди 1 
частоти) з його спектра. Аналзований у статт! метод е альтернативою застосування вконних функщй, що 
пригн!чують, але не виключають ефект «витоку» спектрально! потужност!. На вдм!ну в1д застосування 
вконних функщй, представлений у статт! метод дозволяе практично повнстю вдновити 1стинн! значення 
амплтуди 1 частоти сигналу за даними дискретного спектра. У статт! аналзуеться робота пропонованого 
чисельного методу усунення ефекту «витоку» спектрально! потужност! при разному ступен! шуму вихдного 
перюдичного сигналу. В якост! прикладу взято аналзатор спектра в1броакустичного сигналу турбонаддуву 
судового дизеля, побудований на баз! контролера АТхтега256АЗ. 

Ключов! слова: спектральний аналпз, ефект «витоку» ДПФ, судовий дизель, турбонаддув, 
длагностика. 


Введение 


В сборнике трудов «Искусственный интеллект» за 2012 г. была опубликована 
статья [1], в которой показана возможность определения частотных характеристик 
современного комбинированного судового дизеля путем спектрального анализа вибро- 
акустических полей. В результате применения такого метода можно с высокой точно- 
стью определять частоту вращения ротора газотурбонагнетателя (ГТН) и коленчатого 
вала, а также производить диагностику технического состояния ГТН и оценку нерав- 
номерности цилиндровых мощностей. В качестве анализируемых величин рассматри- 
ваются частоты и амплитуды отдельных гармоник. В этом случае эффект «утечки 
мощности» спектра, построенного в результате ДИФ, приводит к существенным иска- 
жениям в определении частоты, фазы и амплитуды гармоник. 

Практическому решению этой проблемы посвящено не очень много работ и основ- 
ным методом, применяемым на практике в большинстве случаев, является метод оконных 
преобразований исходного сигнала [2]. Проблема в том, что метод оконных преобразова- 
ний не исключает, а лишь уменьшает действие эффекта «утечки» мощности и искажения 
истинных величин амплитуды, частоты и фазы анализируемого сигнала. Другими 
словами, расчетные параметры исходного сигнала, которые необходимо определить 
путем его спектрального анализа, зависят от соотношения частоты этого сигнала 
и частоты дискретизации. Частота дискретизации может быть задана постоянной 
(как принято в большинстве технических систем), а частота анализируемого сигнала 
меняется произвольно. Таким образом, погрешность оценки параметров измеряемого 
периодического сигнала будет не постоянной и меняться вплоть до 50%. Очевидно, 
что без решения этой проблемы ни о каком численном количественном анализе внеш- 
них периодических сигналов с помощью ДПФ не может идти речь. Предложенный в [1] 
численный итерационный метод решает данную задачу. Возникает вопрос, как будет 
работать этот метод в случае наличия шумов в исходном сигнале? 


Определение частотных характеристик и диагностика 
газотурбонагнетателя судового дизеля 


Модуль, определяющий частотные характеристики газотурбонагнетателя судового 
комбинированного дизеля, построен на базе микроконтроллера АТхтега256АЗ (табл. 1). 
Основными факторами, определяющими алгоритмы работы модуля, являются объем 
его оперативной памяти и параметры АЦП. Частота преобразования АЦП достаточно 
велика: 2МУРУ, что на два порядка превышает максимально возможные частоты исход- 
ных сигналов. Объем оперативной памяти, равный 16КВ, ограничивает длину ряда 
исходной выборки для процедуры ЕЕТ - 2048 точек. 
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Таблица 1 — Характеристики микроконтроллера АТхтега256АЗ 


РазН (КБу*е5): 256 КБуе5 Сатега № рисоРо\мег: Уе5 
ГиеГасе: 
Рш СоипЕ: 64 АРС сПаппе[5: 16 Тетр. Капде -40 о 85 
(аед С): 
Мах. ОрегаЧ пд 32 МНЕ АРС Кезо!ийоп 12 О 5ирр!у Саз$: 1.6 © 3.6 
Егеачепсу: (5165): 
СРИ: 8-5 А\УВ АОС 5рееа 2000 Орега#тад \Мо{аде 1.6 © 3.6 
(К$р$): (\Мсс): 
# оЕ ТоисН 16 Апа!од 4 ЕРИ: № 
СВаппе[5: Сотрагафогс: 
Нага\м/аге № Вез15Иуе Тоисй № МРЦ / ММЦ: по / по 
ОТочсН 5сгееп: 
Асаш$Шоп: 
Мах 1/О Рт$: 50 РАС СВаппе[5: 2 Титег$: 7 
ЕхЕ Тпеггир{$: 50 РАС Везо!иНоп 12 ОшриёЕ Сотраге 22 
(55): сПпаппе[5: 
У5В 5рееа: № Тетр. 5епзог: Уе5 ТприЕ Саркиге 22 
СПаппе!5: 
Ц5В ПщеНасе: № Сгурто Епдте: АЕЗ/РЕЗ Р\/М Сваппе[ $: 22 
РГ: 10 ЗВАМ (КБуЕе5): 16 З2КН2 ВТС: Уе5 
тмл (2С): 2 ЕЕРКОМ (Ву*е$): 4096 СаНЬгаеа ВС Уе5 
ОзсШаог: 
АВТ: 7 беЁ Ргодгат \ЕЗ \/аспаод: Уе5 
Метогу: 
СгарЫ!с ЕСО: № Ежегпа! Ви$ 0 
ГеГасе: 
\!аео Ресодег: № ОКВАМ Метогу: № 


В качестве примера объекта диагностирования возьмем ГТН УТК-350 главного 
дизеля 6Г8ОМСЕ. Спектр виброакустического сигнала компрессора УТВ-350 представ- 
лен на рис. 1, на котором показан диапазон интересующих информативных частот: 
от 100 Гц (слева от 1-й гармоники) до 5 КГц (справа от 2-й гармоники). 


147.4 На 
2948 Н2 120 *`60 = 884 трт 


30 50 100 200 300 500 1000 2к ЗК 5К 10000 


Рисунок 1 — Спектр виброакустического сигнала ГТН УТВ-350 


Если производить набор данных с частотой АЦП, равной 50 КГц, то получим 
в выборке из 2048 точек: 
пл! = [2048/50000] х100 = 4.096 периодов. 
112 = [2048/50000] х5000 = 204.8 периодов. 
Соответственно, записанных точек на один период: 
р! = 2048 / 4.096 = 500; 
р2 = 2048 / 204.8 = 10. 
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Это количество периодов и точек на период показано для максимального и ми- 
нимального по частоте исходного сигнала. На номинальной нагрузке частота измеря- 
емого сигнала будет близка (1) к 150 Гц (=10%) и (2) кЗ кГц (+10%). В этом случае 
количество периодов в выборке с частотой 50 КГц будет: 

т! = [2048/50000] х150 = 6.144 периодов. 
т? = [2048/50000] х3000 = 122.88 периодов. 
Записанных точек на один период: 
р: = 2048 / 6.144 = 333.(3); 
р2 = 2048 / 122.8 = 16.68. 

Такого количества периодов измеряемого сигнала в выборке и точек на один период 
вполне достаточно для корректного отображения информативных гармоник в дискрет- 
ном спектре, построенном с помощью БИФ. Во всех случаях гарантированно выпол- 
няется требование теоремы Котельникова [2], [3]. 

Информативные гармоники (1 и 2, рис. 1) в спектре [4] показывают следующее: 

— лопаточную частоту компрессора ГТН УТВ-350 (2), определенную с помощью 
ограничений [1]. Ее деление на количество лопаток позволяет с большой точностью 
непрерывно определять текущую частоту вращения ротора ГТН, которая находится 
в прямо пропорциональной зависимости от эффективной мощности двигателя: 

ПТОК = 2948 Гц /20 * 60 = 8844 грт; 

— центральную гармонику (1). Ее амплитуда показывает колебания ротора на основ- 
ной частоте. Повышение амплитуды этой гармоники является признаком нарастающего 
динамического разбаланса ротора ГТН и, соответственно, комплексным диагностиче- 
ским признаком. Динамическая балансировка ротора ГТН может нарушаться в случаях 
дефектов подшипников, либо в случае засорения соплового аппарата, лопаток и проточ- 
ной части ГТН: 

УГОК = 2948 Гц / 20 = 147.4 Гц. 


Устранение эффекта «утечки» в дискретном спектре 


Эффект «утечки» является следствием конечности анализируемой временной реа- 
лизации и ее дискретного представления через АЦП. Он считается одной из главных 
погрешностей дискретного преобразования Фурье. 

Если частоту сигнала представить в виде 

у = М/Т, 

где Т — период сигнала, а М =п + с, где п — целое, а 0 <с < 1, то максимальные 
искажения амплитуды, частоты и фазы центральной гармоники и утечка мощности в со- 
седние будет наблюдаться при о = 0.5 [1], [2]. 

Таким образом, если анализировать параметры стационарного в течение выборки 
сигнала по его спектру, т.е. по центральной гармонике, то полученные амплитуда, 
частота и фаза будут искажены в случае нецелого количества отсчетов сигнала на его пе- 
риод. На практике для дискретной записи сигналов используют АЦП с выбранной 
и зафиксированной частотой дискретизации. Понятно, что число отсчетов на период 
никогда не будет целым, а значение с будет меняться от 0 до 1 в зависимости от соб- 
ственной частоты измеряемого сигнала, при этом будет меняться точность оценки па- 
раметров сигнала по центральной гармонике. 

Для уменьшения эффекта «утечки» наиболее широкое распространение получил 
метод оконных преобразований: для уменьшения разрывов на краях ряда с целью ослаб- 
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ления утечки нужно уменьшить амплитуду сигнала возле краев. Такое масштабирование 
осуществляется в ходе умножения реализации на окно специальной формы х” =х, .И/’() 
(окно Хемминга, Ханнинга, Блекмена, Кайзера, и т.п.). В результате применения окон- 


ных функций спектр исходного сигнала может сильно исказиться, но при этом умень- 
шится зависимость амплитуды, частоты и фазы от значения о. Это значит, что централь- 


ную гармонику х” можно использовать для приблизительной оценки параметров сигнала 


с некоторой постоянной погрешностью, которую можно учесть. Точнее можно устра- 
нить эффект «утечки» численным методом, основанным на обработке комплексных 
результатов ДПФ. В [6] предложено уточнить частоту т, фазу ф, и амплитуду А исход- 
ного сигнала по значениям двух максимальных гармоник в спектре Х,,Х,,,. Для этого 
предлагается численно решить систему комплексных уравнений: 


|2 (т,Ф), /2(т,Ф) ии = |х, р 


Аз (Е (т,ф),) = Ат8 (Ху) 
где параметры К-й гармоники: 


(0 


Х, =Ке, + 7 1щ, ; 
Я, = Ма", 


А, = уе Пить , 


Тп 
ф, зат 


)= 4"8(Х,). 


ый 
Коэффициенты гармоник можно представить в виде: Х, = (А, /2)7(т,Ф),, 
где 2(т‚ф), — комплексная функция, не зависящая от амплитуды, но зависящая от ча- 


стоты и фазы: 


— р29ат-®) _ 1 , ото ы] 
— > —/Ф 
Е(т,ф), =е 217 (т—®) +е 2 (тк) 
е м -| е № -1 


Систему уравнений (1) решается в том случае, если гармоники слева и справа 
от центральной не равны нулю, т.е. больше заданной малой величины д: 
Роль. 
Система (1) решается методом итераций. После завершения итераций амплитуда 
сигнала определяется из соотношения: 
А, = А65(2Х,/7(т,ф),). 


Восстановление параметров зашумленного сигнала 


На рис. 2 показано решение системы (1) для чистого сигнала и сигнала с добав- 
лением 5% и 10% шумов стандартной функцией Вапдот(). Для всех исследованных 
случаев требовалось не более 6-и полных итераций для обеспечения заданной погреш- 
ности менее 0.5% по частоте и фазе. Сходимость обеспечивалась во всех случаях, даже 
при уровне шумов 20 %. Причем, среднее количество итераций -— 5 — наблюдалось 
для любых значений о, от минимально выраженного эффекта «утечки» при с - 0, 
до максимального при с - 0.5. Амплитуда и частота во всех случаях были восстанов- 
лены до заданного в исходном сигнале значения (А1=0.8, \1=5.0) с точностью до 5-0 
знака после запятой. 
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* зе б2ти/Т + 2 145/484 > 
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Рисунок 2 — Окно программы демонстрации ДПФ и решения проблемы «утечки» 


При решении системы (1) используется много повторяющихся участков, которые 
удобно выделить в подпрограммы. Решение системы (1) не связано с выделением 
дополнительной памяти для хранения объемных массивов данных и расчетных коэф- 
фициентов, как в случае с БПФ. В связи с этим алгоритм может быть запрограмми- 
рован на современном О$Р контроллере, реализующем БПФ. 


Заключение 


Несмотря на итерационное численное решение системы (1), эта процедура восста- 
новления дискретного спектра незначительно увеличит общее время расчета и позволит 
получить не только спектр сигнала, но и восстановленное значение основной частоты 
и амплитуды исследуемого сигнала. Виброакустический метод определения частотных 
характеристик надувного судового двигателя имеет ряд преимуществ перед традицион- 
ными. В первую очередь это стоимость и простой способ инсталляции датчиков. Кроме 
того, постоянный анализ спектра виброскорости/ускорения блока цилиндров двига- 
теля и ГТН дает возможность анализировать техническое состояние двигателя и газо- 
турбонагнетателя. 

Как выяснилось, предлагаемый метод восстановления дискретного спектра ока- 
зался устойчив к шумам в исходном сигнале. Это дает возможность его широкого приме- 
нения в составе систем непрерывного мониторинга и диагностики современных судовых 
дизелей и их систем газотурбонаддува. 
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